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In the dilute limit, the properties of fermionic lattice models with short-range attractive interac- 
tions converge to those of a dilute Fermi gas in continuum space. We investigate this connection 
using mean-field and we show that the existence of a finite lattice spacing has consequences down to 
very small densities. In particular we show that the reduced translational invariance associated to 
the lattice periodicity has a pivotal role in the finite-density corrections to the universal zero-density 
limit. For a parabolic dispersion with a sharp cut-off, we provide an analytical expression for the 
leading-order corrections in the whole BCS-BEC crossover. These corrections, which stem only from 
the unavoidable cut-off, contribute to the leading-order corrections to the relevant observables. In 
a generic lattice we find a universal power-law behavior n 1 ^ 3 which leads to significant corrections 
already for small densities. Our results pose strong constraints on lattice extrapolations of dilute 
Fermi gas properties. 
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I. INTRODUCTION 

One of the most attractive features of ultracold Fermi 
gases is the universal physics arising from their intrin- 
sic dilution. A Fermi gas of volume density p is dilute 
when kpRo <C 1, where kp = (37r 2 p) 1 / 3 is the Fermi mo- 
mentum and Rq is the range of the interaction. When 
this condition is realized the interactions are controlled 
only by the s-wave scattering length a s and the details 
of the potential become irrelevant . Therefore the physics 
is controlled by the inverse gas parameter r\ = (kpa s )~ 1 
and a universal phase diagram can be drawn, regard- 
less of the nature of the fermions involved [l[. More- 
over, Fano-Feshbach resonances @ allow to tune the in- 
teraction from weak- to strong-coupling realizing (i) The 
crossover from the Bardecn-Cooper-Schricffcr (BCS) su- 
perfluidity (?/ — > — oo) to a Bose-Einstein Condensation 
(BEC) of preformed pairs (t) -> +oo)g| and (ii) The uni- 
tary limit in which the scattering length diverges, and 
the system has no characteristic length scale other than 
kp 1 leading to an extra universality which attracted a 
remarkable interest 0]. 

In a different framework, the physics of superconduc- 
tors beyond the BCS regime, including the high-Tc ma- 
terials, has triggered the investigation of lattice models 
with attractive interactions. The typical electronic densi- 
ties in these materials are far from the dilute limit, being 
instead closer to the so-called half-filling regime in which 
the number of fermions equals the number of lattice sites. 
However a BCS-BEC crossover takes place also in this 
regime [5| even if no universality is expected. 

In this work we focus on the connection between these 
two different realizations of the BCS-BEC crossover. The 
main effect of the lattice is to lift the Galilean trans- 
lational invariance introducing a characteristic length 
scale, the distance between two lattice sites, or lattice 
spacing. The reduced translational invariance makes the 



fermionic spectrum bounded from above or, in other 
words, introduces an ultraviolet energy cut-off. 

When the number of fermions is of the same order of 
the number of sites the underlying lattice influences the 
properties of the fermionic gas substantially. Reducing 
the density, the lattice becomes less and less relevant, 
and in the dilute limit the properties of lattice fermions 
are expected to converge to the universal features of di- 
lute gases. The evolution of the properties of a lattice 
model by reducing the density is an interesting topic per 
se, also in light of the possibility to study the BCS-BEC 
crossover in optical lattices at various particle densities 
@, but it can also be used to benchmark and test theo- 
retical studies in which the properties of dilute gases are 
obtained through a low-density extrapolation of lattice 
models, as done in recent lattice Quantum Monte Carlo 
(QMC) approaches to the unitary Fermi gas 

In a previous work [l8l ]. we used static and dynamical 
mean-field theory [1^ at the unitary limit. We proved 
that finite-density corrections to the dilute unitary limit 
can be large even at small densities and that the depen- 
dence on the specific lattice is hardly washed out at any 
finite density. Therefore a safe extrapolation to the dilute 
limit is in general extremely difficult. 

In this work we extend the analysis in the whole BCS- 
BEC crossover and we focus on the origin of the finite- 
density corrections. In order to gain some analytical un- 
derstanding, we use a static mean-field theory, which al- 
lows to treat on the same footing the lattice and contin- 
uum cases. The simplicity of the approach helps us to 
identify the origin of the leading density-corrections to 
the dilute limit together with their functional dependence 
on the inverse gas parameter r\ in the whole crossover. 

The paper is organized as follows: In Sec. |IIJ we define 
the connection between lattice and continuum models in 
terms of the relevant variables in the dilute regime. Then 
we introduce in Sec. IIIII the mean-field approach and 
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we recast it in terms of these variables. Results for the 
unitary system and in the whole crossover are shown in 
Sec. IIVI Finally we draw the concluding remarks in Sec. 
IVIl together with the consequences of our findings over 
lattice extrapolation of the BCS-BEC crossover in dilute 
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TABLE I: DOS summary (H=l = m = V) 



II. THE LATTICE PATH TO THE DILUTE 
LIMIT 

In this section we establish formally the connection be- 
tween lattice models and dilute gases in continuum space. 
First we identify the lattice counterpart of the inverse 
gas parameter 77, then we discuss the role of the specific 
dispersion of the lattice and of the ultraviolet cut-off. Fi- 
nally we cast some general considerations in the form of 
a phase diagram in the relevant parameter space. 



A. Interaction: the scattering length 

As briefly discussed in the Introduction, for a dilute 
Fermi gas interacting via a short-range potential the 
only relevant parameter to describe the interaction is 
the s-wave scattering length a s and the properties of 
the system arc defined by the inverse gas parameter 
i] = (fcpa s ) . In this subsection we define the lattice 
versions of a s and rj. 

We model the lattice system with an attractive Hub- 
bard model 



H = -t 
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(1) 



where t is the hopping parameter, U < is the onsite 
attraction and /1 is the chemical potential which fixes the 
density of particles per lattice site n in a grandcanonical 
ensemble. 

Since the scattering length a s is a two-body property, 
it can be evaluated from the sole knowledge of the in- 
terparticle potential. a s is proportional to the limit of 
vanishing energy and quasi-momentum of the the two- 

i.e. 
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where m is the bare effective mass of the fermions in the 
lattice at hand and I is the lattice spacing. 

For a Hubbard model T l £\ dy (u , q) = [U' 1 - 
H latt (w, <f)] _1 , where H latt is the polarization bubble 
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the momentum integral extends over the first Brillouin 
zone and the zero energy is chosen to have non-negative 
energies and eg = 0. 
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where U c is defined by [II Jat *(0, 0)] 1 = U c and, using 
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In the expression above we emphasized that the lattice 
density of states (DOS) Di att is finite only below a cut-off 
A. The scattering length therefore diverges at U c which 
defines the unitary limit on the lattice. It is worth to no- 
tice that for \U\ — > 00 a s converges to a finite value of the 
order of the lattice spacing /. Hence even for divergent 
interactions the finite lattice spacing implies a non-zero 
size of the Cooper pairs. 

In a lattice model it is natural to describe the system 
is terms of the number of particles per lattice site n (also 
called filling factor), n/l 3 is the number of particles per 
unit volume and we can define the Fermi energy Ep and 
the Fermi momentum ^ of a non- interacting homoge- 
neous gas with the same density and mass m as 
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(01 and © define the lattice analogue of the inverse gas 
parameter 77 = (fcjro, ) — 1 . In the following we will study 
the attractive Hubbard model along lines with constant 
77 which span the whole BCS-BEC crossover. 



B. Lattice dispersion and the role of the cut-off 

As outlined in the introduction, an unavoidable con- 
sequence of the broken Galilean invariance is the exis- 
tence, in any finite dimension, of an ultraviolet cut-off A 
in the dispersion. This corresponds to the finite band- 
width of lattice systems, as opposed to the unbounded 
parabolic dispersion of particles in free space. We will 
show that, even if the cut-off occurs at energies much 
larger than the Fermi energies corresponding to low- 
density (Ep oc n 2 / 3 ), it contributes to the leading-order 
corrections in the density to the universal dilute limit. 
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In order to properly reproduce the dilute limit of 
fcrmion gas with mass m in three-dimensional homo- 
geneous space we have to consider a class of lattice 
dispersions which share the same low-energy behavior 



of the free DOS D free (e) 
Aatt(e) oc y/e (or e k oc k 2 ). 
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In order to identify how the lattice and continuum 
physics connect, we study the following three models 
(shown in Fig. [I}, all sharing the correct low-energy be- 
havior, but substantially different at high-energy 



A semicircular DOS 
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where D is the half-bandwidth and A = 2D is the 
cut-off and D = (4tt) 2 / 3 -C. 

• A first lattice gas model (LG\) with a sharp cut-off, 
namely a system in which the DOS is that of free 
fermions up to a sharp energy cut-off A^Gi which 
enforces the presence of the lattice and ensures the 
correct normalization of the DOS per lattice site. 
Above the cut-off the DOS is identically zero. 

• A second lattice gas model (LG2) with a smoother 
cut-off. The dispersion coincides again with that of 
free fermions, but the lattice periodicity is enforced 
restricting the momenta to the first Brillouin zone 
of a cubic lattice (k € [ — tt/Z, 7t/Z) 3 ). In this way the 
DOS coincides with the free one only up to AiG 2 /3, 
then smoothly decreases and vanishes at Alg 2 - 

By comparing these three models, we can clarify the 
roles of the cut-off A and the effect of the different be- 
havior of the DOS'S at the various energies. 

In order to properly compare the different lattices, we 
evaluated according to Eq. §5§ the critical value U c for 
each lattice. Results are summarized in Tab. HI together 
with the value of the cut-off and other quantities defined 
in the following. 

We do not study the cubic lattice, which has been 
shown in Rcf. [lj§] to introduce larger finite-density ef- 
fects than the models defined above. 



C. Approach to the dilute limit 

Once the mapping (U,n) — > (a s ,kp) is established via 
Eqs. (0| and ([5]), we can consider the parameter space of 
the lattice model at T = in terms of the inverse gas pa- 
rameter 77, which controls the physics in the dilute limit. 
By expressing 77 as a function of the coupling U and of 
the density n and inverting this relation, we can calcu- 
late the trajectories for fixed 77 in the (|J7|,n) parameter 
plane of the attractive Hubbard model, i.e. 
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FIG. 1: (Color Online) Comparison between the different lat- 
tice DOS under investigation: Semicircular (full red line), 
Lattice Gas 1 (dashed blue line), and Lattice Gas 2 (dashed- 
dotted green line). For reference the DOS per lattice site 
l 3 Df ree for free fermions in three-dimensional space is shown 
as a dotted black line. We set h = I = m = 1. 



where a 



is a dimcnsionless constant 



depending on the lattice at hand. In Fig. [2] we plotted 
these trajectories for several values of the inverse gas pa- 
rameter 77 in the BCS-BEC crossover for a semicircular 
DOS. These trajectories give a pictorial representation 
of how the physics of a lattice model at finite density 
evolves towards a continuum system in the dilute limit. 

It should be clear, however, that this diagram is not 
meant to provide information about the density where 
we expect to get rid of non-universal lattice features and 
access the universal properties of dilute gases. This is 
particularly evident for the trajectory at unitarity 7/ = 0, 
where, even though the correct value of the interaction is 
U c for any density, the universal properties of the unitary 
Fermi gas can be recovered only at small density. 
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FIG. 2: (Color Online) Trajectories with constant 77 in 
the (|[/|,n) plane of an attractive Hubbard model for 77 = 
-10,-1, -0.1 (BCS side), r? = (Unitarity) and 77 = 0.1, 1, 10 
(BEC side). The results shown concern the semicircular DOS 
(see Sec. [H]) but a similar picture is expected in the other 
cases. 
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As shown in Fig. $Z§, the trajectories converge at low- 
density towards the critical value U c . While at large den- 
sities the crossover is wide in U, in the dilute region ev- 
erything happens in a very narrow region around U c , and 
the values far from U c are safely in the BCS and BEC 
regimes, respectively. 

We expect that the convergence to the dilute regime 
is slower on the BEC side, where the property of the su- 
pcrfluid state are more sensitive to the presence of the 
lattice. Indeed, in the strong coupling regime for the 
lattice model ^ 3> 1 the size of the superfluid pairs is 
bounded from below by the lattice spacing. In this regime 
the fermions form strongly bound pairs with a large ef- 
fective mass Mb oc ^ 3> m (heavy bosons region in 
Fig. [2]). Despite of the hard-core constraint arising from 
the Pauli principle being irrelevant at low-density, even 
a single pair can only move through virtual ionization 
due to the presence of the lattice. This implies that the 
reduced translational symmetry of the lattice is relevant 
at any finite density in this regime. 

It is useful to notice that, due to the fact that the scat- 
tering length cannot be arbitrarily reduced in the lattice, 
the trajectories for 77 > have an horizontal asymptote 
for large \U\, implying that large positive values of rj can- 
not be accessed in the lattice above a given density. 

On the other hand, the BCS regime is not expected to 
be qualitatively dependent on the density allowing for a 
smoother extrapolation to the dilute limit. This expec- 
tation is confirmed by the results that we present in the 
Section [TV] 



III. MEAN-FIELD APPROACH 

After the general considerations above, we analyze 
quantitatively the convergence of the lattice models to a 
dilute Fermi gas in continuum space by means of a static 
mean-field approach. As shown by Leggctt [20] for the di- 
lute gas and by Nozieres and Schmitt-Rink for the lattice 
model [2l| , the mean- field approach provides a reasonable 
description of the ground state properties through the 
whole BCS-BEC crossover, being exact both in the BCS 
and BEC limit and interpolating between them at inter- 
mediate coupling. Within this simple approach we can 
deal with continuum and lattice system on equal footing 
without introducing any bias. Another advantage over 
more accurate methods is that it can be implemented 
directly in the thermodynamic limit, getting rid of finite- 
size effects, and in the grand canonical ensemble, which 
allows to address arbitrarily small densities. 

We implement a mean-field approach for the attractive 
Hubbard model and we follow trajectories of fixed ?y in 
the parameter space. For vanishing density we expect 
to recover the results of Leggett's mean-field theory for 
dilute Fermi gases [20j | . 

The static mean-field equations for the Hubbard model 
in the superfluid phase are easily obtained by consider- 
ing the Hamiltonian (TTJ) and decoupling the interaction 



term both in the particle-hole and in the particle-particle 
channel and introducing a gap parameter A = U(ci^Ci y \) . 
The self-consistency equations read 
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where ((e) = (e — //), 7(e) = y/(e — /1') 2 + A 2 and the 
Hartree term is included in the shifted chemical potential 
p! = fi — (Un/2). The ground state energy per lattice site 
E is given by 
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where Eq = Un 2 /4 + A 2 /U. The self-consistency equa- 
tions © and (|TU)) have been already studied in [21j as a 
function of the interaction U at fixed density. The focus 
of this work is, however, the evolution towards low den- 
sity and the connection for n — > to the Leggett mean- 
field approach in the continuum [20j . In this way we can 
use them to gain insight on the finite-density corrections 
to the dilute limit in a lattice. 

By solving numerically the mean-field equations for the 
different lattice models under consideration we can assess 
that our method reproduces Leggett's results as derived 
in Ref. [22| in the whole BCS-BEC crossover and, most 
importantly, we quantify the deviations from the dilute 
limit as a function of the density and interaction, together 
with their dependence on the specific lattice at hand. 

By inverting the expressions of the Fermi energy Ep in 
Eq. ^ and of the inverse gas parameter rj(U, n) in Eq. 
©; wc can recast Eqs. (|9I10[) in terms of the relevant 
variables for the dilute limit, obtaining an equivalent set 
of equations for the renormalized gap parameter A 
and the renormalized chemical potential jl 

After some algebra one obtains 
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where x = e/Ep, D(x) 



~ _ nU v (n) 
t L 1Ep 



Diatt(E F x) ~, 
I 3 D j ree (E F ) ' A 4 

is the renormalized chemical potential including the 
Hartree contribution, A is the energy cut-off which de- 
rives from the broken translational invariance and corre- 
sponds to the lattice bandwidth in the solid state frame- 
work. We remark that we are not assuming that the den- 
sity n is small at this step, but we just reformulated the 
mean-field equations in a lattice in terms of the relevant 
variables in the dilute limit. 

It is easy to verify that for n — > 0, i.e. Ep — > 0, and 
fixed 77, our mean-field equations reduce asymptotically 
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to Leggett mean- field equations j2Cj. In this limit both 
the gap and the chemical potential are proportional to 
Ep through coefficients which are only function of rj. On 
the other hand, as we shall discuss in the following, for 
any finite density fl and A retain an additional depen- 
dence on the density according to Eqs. (|12ll3p . 

The corrections to Leggett theory arise from two dif- 
ferent sources: (i) the different DOS of the lattice mod- 
els and (ii) The Hartree term, which does not appear in 
Leggett approach. 

As far as the DOS is concerned, we can disentangle 
two effects, namely the simple fact that at any energy the 
lattice DOS can be different from the continuum coun- 
terpart, and the existence of the cut-off energy A which 
directly stems by the existence of a finite lattice spacing. 
Even though it is possible to choose the lattice DOS to be 
identical to the target DOS at low-energy, as in the case 
of the lattice gas models introduced in Sec. Ill BL there is 
no way to avoid the breakdown of Galilean invariance. 

The Hartree term, absent in the Leggett approach, 
vanishes as n 1 ^ 3 once rescaled with Ep. Even if it can be 
easily reabsorbed in a shift of the chemical potential, it 
plays a leading role in the asymptotic corrections to the 
dilute limit, as we will clarify in the next section . 

Finally we notice that, keeping the Fermi energy Ep 
fixed (i.e. for fixed density n), the same asymptotic equa- 
tions are achieved if the cut-off A is divergent, corre- 
sponding to a vanishing lattice spacing. This alternative 
point of view is the analogous of taking the zero-range 
limit of a continuum model instead of reducing the den- 
sity of a system with a finite range b. 

IV. RESULTS 
A. Unitary Limit 

As briefly mentioned in the Introduction, dilute Fermi 
gases evolve continuously as the interaction is tuned from 
the BCS limit (77 -> -oo) to the BEC limit (77 -> 00) 
and a crossover takes place in the region around 77 = 
0. We start our analysis of the mean-field results from 
the unitary limit, where the scattering length diverges, 
corresponding to 77 = 0, or cquivalently U = U c . 

This limit has attracted an enormous interest because 
the divergence of a s makes any length scale other than 
kp 1 irrelevant leading to an extra universality. Precisely 
for this reason, most of studies attempting to extrapolate 
the properties of dilute gases from a lattice perspective 
have concentrated in this limit. 

For the sake of definiteness we start with the results for 
the semicircular DOS, which presents some advantage in 
the numerical solutions. In Fig. [3] we show the rescaled 
gap parameter A, chemical potential fl and internal en- 
ergy £ = 3/5n_E F as a function of the lattice density n. All 
the quantities correctly converge to Leggett's mean-field 
values for the dilute gas (A^ ~ 0.68 and /it, = £l — 0.59) 
as n — >• 00. Two important facts emerge: All the curves 



approach their limiting values from below with a diver- 
gent slope, and the finite-density corrections are hardly 
negligible unless the density is extremely small. For ex- 
ample, in order to have a deviation smaller than 1% for 
the three observables, we need n < 10~ 3 . 

A free fit of the data in the low-density region (n < 
10 -3 ) shows that the leading asymptotic corrections to 
all the observables under consideration are proportional 
to n 1 / 3 for any lattice studied and in the whole BCS-BEC 
crossover (i.e. for any 77). The same behavior has been 
reported for the critical temperature in Refs. [Iol - [l3j |. In 
TableUwe report the coefficients of the n 1//3 for the three 
observables a M , oa, and a^. 

The observed scaling is consistent with the hypothe- 
sis that the rescaled observables are analytic functions of 
kpl for every value of 77. The linear term in 77 clearly re- 
sults in a n 1 / 3 correction at small densities. In Fig. 0]we 
plot the same quantities as a function of kpl for the three 
lattice models defined above. Even if the corrections to 
the universal dilute limit are of the same order of magni- 
tude for any lattice, the semicircular DOS gives smaller 
corrections in fl and £ within our mean-field approach. 

In the following we discuss the different contributions 
to these leading order corrections and the actual role of 
the properties of each lattice. 

The simplest term to analyze is the Hartree correc- 
tion to the rescaled chemical potential, which at unitar- 
ity (U = U c ) is simply 8(1 = |^ and contributes to the 
leading-order corrections proportional to n 1 / 3 . The nu- 
merical value is 5fx = jjjP; = — 0.171(kpl) for the LGi 
and it is of the same order for all the lattices. On the 
other hand the calculated coefficient a M is negative, but 
significantly smaller in magnitude (see Table HJ), suggest- 
ing that the large negative contribution from the Hartree 
term has to be compensated by a positive contribution 
which necessarily arises from the difference between the 




FIG. 3: (Color Online) Rescaled chemical potential jl, gap A, 
and ground state energy £ for the semicircular DOS case along 
the unitary trajectory (77 = 0) as a function of the density n. 
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FIG. 4: (Color Online) Rescaled chemical potential jl, gap A 
and ground state energy £ at unitarity (77 = 0) as a function 
of the lattice dilution parameter k F l for the semicircular DOS 
(red circle), the lattice gas 1 (blue squares), and lattice gas 2 
case. 



In the following we discuss the origin of these lattice- 
dependent positive contribution, disentangling a contri- 
bution arising from the presence of the cut-off A and a 
more subtle contribution due to the fact that the lattice 
DOS may be different from the continuum one at finite 
energy. 

In order to disentangle these contributions, we consider 
the LGi model, in which the DOS is identical to the free 
one for all energies below the cut-off, which is obviously 
the only source of corrections. 

We can rewrite the mean-field equations p"2")) and (fT3")) 
for the LG\ model in a compact form as 4/3 — A A (jl' , A) 
and 77 = B A ([i' , A), where we introduced 



A A (fl,A) = / dxy/x 
Jo 
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x — fl 
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(14) 



(15) 



and A = ^ . For n — > A diverges and Leggett's equa- 
tions are recovered. By expanding the equations above 
for large values of A and linearizing around the Leggett 
solution (fJ-L, A-l), we can derive an analytic expression 
for the leading-order lattice corrections, i.e., 
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where a = E F /n 2 ^ 3 and C\ and C2 are constants calcu- 
lated in the zero-density limit which do not depend on 
the lattice 
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and Det[J] F is the determinant of the Jacobian matrix 
J = g^^j evaluated for Ji = Jil and A = A^. At uni- 
tarity Ci ~ 0.4170 and C 2 ~ -0.0775. Eqs. (Hlj) and 
(|17[) are exact expressions for the leading-order density 
corrections arising from the very presence of a finite cut- 
off energy in the LG\ case. In the case of the chemical 
potential this is the term that partially compensates the 
Hartree shift, while it gives the full leading-order correc- 
tion to the order parameter. 

These results unambiguously prove that the presence of 
a cut-off determines n 1 / 3 corrections to the observables, 
and there are no self-evident choices of the lattice that 
can cancel these corrections that ultimately arise from 
the broken Galilean invariance. 

Obviously for the other lattices, other corrections due 
to differences between the free DOS and the specific lat- 
tice DOS are expected. Indeed the expressions ([T^]) and 
([XT]) (evaluated introducing the proper values of U c and 
A shown in Table H]) do not reproduce the calculated co- 
efficients of the n 1 / 3 corrections for the semicircular DOS 
and LG\ models. This implies that also the behavior of 
the DOS at finite energy contributes to the same order 
of the cut-off (and of the Hartree shift in the case of the 
chemical potential) and therefore leads to further n 1 / 3 
corrections. 

Werner and Castin [23[ have shown for a gas in the con- 
tinuum space that the leading-order density to the dilute 
limit assume a universal behavior in terms of the effec- 
tive dilution parameter kpR e , where the effective range 
R e is defined through the low-momentum k expansion of 
the scattering amplitude / [24[, i.e. 
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ik~R e k 2 /2 
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One could expect a similar universal behavior on a 
lattice if the asymptotic corrections to the dilute limit 
were determined by the low-energy scattering proper- 
ties alone. In Fig. [5] we plot the usual observables as 
a function of k F R e using the calculated values of R e for 
our lattices. Our results do not show any universality, 
confirming that the broken Galilean invariance, i.e. the 
simple existence of a cut-off in the dispersion, plays a 
key role in the leading-order corrections, which are not 
determined only by the low-energy scattering properties. 
In particular within MF the leading order corrections are 
of order n 1 / 3 also for the lattice introduced in Ref. [TtJ 
and designed in order to have R e = (not shown). 
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f<] 0.65 




FIG. 5: (Color Online) Rescaled chemical potential jl, gap A 
and ground state energy £ at unitarity (77 = 0) as a function 
of the lattice effective dilution parameter kpRe for the semi- 
circular DOS (red circle), the lattice gas 1 (blue squares), and 
lattice gas 2 case. 



Our results also show that there is no easy way to 
eliminate or systematically the non-universal corrections 
to the dilute gas behavior even at very small density. The 
choice of the lattice can be important only to quantita- 
tively reduce the finite-density correction, but no lattice 
can remove the effects of the reduced translational invari- 
ance. 



B. Crossover Region 

In this section we extend our analysis beyond the uni- 
tary point. Therefore we solved the mean- field equations 
as a function of the density along the trajectories at con- 
stant 77 = in the whole crossover spanning the region 
— 2 < 77 < 2. We start from an overview of the results for 
the SC DOS. 

For every value of 77 we correctly recover Leggett's re- 
sults in the dilute limit. Since the limit values A^ and 
fiL strongly depend on 77 and vary over several orders 
of magnitude in the crossover, it is more instructive to 
plot the relative deviations D (upper panel of Fig. [5]) and 
analogously (lower panel) as a function of the density 
for several values of 77. 

A first general conclusion is that the approach to the 
asymptotic values is slow for all values of 77, confirming 
that, at any value of the effective interaction, a reliable 
extrapolation to the dilute limit requires a reasonably 
small density, whose value depends however on 77. 

As expected from the general considerations in Sec. 
HI! the slowest convergence is observed on the BEC side 
?7 > where the corrections are much larger than for the 
unitary point. The large values of D A observed deep in 



the BCS regime are essentially due to a vanishingly small 
value of the order parameter, and the same apply for 
when 77 = 0.6 since the chemical potential changes sign 
close to this point. 

The gap parameter A is in general less sensitive to 
finite-density effects than the chemical potential. In the 
asymptotic low-density region, the deviation D A (n) is a 
monotonically increasing function of 77 for fixed n. The 
gap parameter is overestimated coming from the BEC 
side of the crossover and underestimated coming from 
the BCS side, even though the gap corrections along the 
unitary trajectory 77 = are more similar to the BCS 
region, i.e., they converge from below. In the same range 
of densities, is negative for every value of 77 we inves- 
tigated, i.e., finite density reduces the rescaled chemical 
potential at any 77. 

We can extend our analysis of the leading-order cor- 
rections to the dilute limit. We start from the Hartree 
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8 



shift of the chemical potential Sft, = 2 Ef wnere we 
emphasized the dependence of \U\ on n along the trajec- 
tories at constant rj. Since U v (n) differs from U c only 
at the order n 1 / 3 (Sec Eq. the leading-order correc- 
tion (ex ) due to this term is identical to the one at 
unitarity for every value of 7). 

This observation implies that the analytic expressions 
for the leading-order lattice corrections due to the cut-off 
in the LG\ model given by Eqs. (|16p and (jTTJ) arc valid in 
the whole crossover, provided that we evaluate C\ and C2 
using /} l (jj) and A p (77) which solve the Leggett equations 
for the given 77. In Fig. we show the results for a^(rj) 
and a&(7]) in the crossover region. 

In the BCS limit 7/ — > — 00 the contribution due to 
the cut-off is vanishingly small, C\ — > + and C2 — > 0~. 
Therefore a M — > —0.171, given by the Hartree shift, and 
aa — > _ . The vanishing of the correction due to the 
cut-off can be understood by considering that in the BCS 
limit pairing only involves an exponentially small region 
of energies around the Fermi energy Ep <C A and the 
cut-off region gives an exponentially small contribution. 

In the case of the LG\ model the Hartree shift and the 
cut-off contribution alone determine the leading-order 
correction, but on a generic lattice, like, e.g., the semi- 
circular DOS model, we can have contributions arising 
from the differences in the functional form of the density 
of states. 

Computing explicitly these corrections for the semi- 
circular DOS in the BCS regime, we find that they are 
of order (kpl) 2 oc n 2 ' 3 . This observation shows that at 
weak-coupling, the difference in the DOS does not con- 
tribute to the leading-order corrections. On the other 
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FIG. 7: (Color Online) Asymptotic slope of the lattice cor- 
rections to the chemical potential a^ Gl and to the gap pa- 
rameter a^ Gl for the LGi model as a function of the inverse 
gap parameter in the BCS-BEC crossover. The inset show a 
magnified view of the the region around the unitary point. 



hand, as shown in our results about the unitary point, 
this is not true outside the BCS regime. 

Turning back to the LG\ model, we find that beyond 
the BCS regime C\ and Ci are always sizeable, generaliz- 
ing the result obtained at unitarity, that the cut-off plays 
a crucial role in determining the leading-order corrections 
even if A 3> Ep. This is essentially due to fact that out- 
side the BCS regime pairing involves a wider region of 
energies whose size increases going to the BEC limit. As 
expected, in the BEC regime the size of the corrections 
explodes, being several order of magnitude larger already 
for moderately large values of 7]. As evident from the in- 
set in Fig. is always negative and has a rather 
broad maximum slightly shifted towards the BEC side. 
aA always increases with 77 being negative in the BCS 
limit and positive in the BEC limit. It is interesting to 
notice that the unitary point T) = docs not show any 
special behavior with respect to the finite-density lattice 
corrections. 



V. OUTLOOK 

Our analysis shows that the presence of a lattice is 
hardly irrelevant even at very small density. Moreover 
we have shown that there is no easy shortcut to over- 
come this poor convergence as a function of density. Even 
what we label lattice gas 1, whose DOS differs from the 
bare one only for the unavoidable presence of a cut-off, 
has sizeable leading-order corrections, and the reduction 
of the effective range can only reduce the quantitative 
effects. 

These considerations are crucial in light of the explo- 
sion of lattice approaches to BCS-BEC crossover in di- 
lute Fermi gases,mainly based on QMC simulations @- 
[Tsj ] and DMFT [25j]. We remark that our analysis only 
refers to lattice simulations and it does not pretend to 
predict the behavior of low-density corrections in contin- 
uum space. 

QMC approaches for attractive models provide very 
accurate results (or at the very least upper bounds) for 
finite-size lattice models at finite density n. Our point is 
that, even for exact methods, a great care is needed to 
extract reliable information on the actual dilute system 
in the thermodynamic limit. 

Even after a careful extrapolation to the thermody- 
namic limit is performed in order to eliminate finite-size 
effects, a further extrapolation of the dilute limit n — > 
is still required to access universal and therefore model- 
independent features in the dilute regime. 

Precise results require therefore simultaneously large 
lattices and large number of particles. We argue that 
one of the reasons why current QMC results are scattered 
around the most precise experimental estimates (see e.g. 
[26l - |29l ] ) rather than converging to them, is that non- 
universal lattice features are still present at the densities 
considered in these approach and not properly taken into 
account. 
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Our analysis clearly shows that choosing a given den- 
sity and assuming that it is small enough for the prop- 
erties of the system to coincide with the dilute limit can 
lead to substantial deviations, especially if the densities 
are of order n ~ 0.1 as in Refs. At these densities 

the deviations from the dilute gas behavior can not be 
considered as small (see e.g. Figs [3] and [6J , especially if 
the aim is the quantitative evaluation of the dilute gas 
properties, like, e.g., the existence of a small pseudogap 
at unitarity. The agreement with recent experimental re- 
sults [26| is most likely resulting from a cancellation of 
residual finite-size and finite-density effects. 

While the need of a careful extrapolation to the di- 
lute limit has been realized by different authors (see e.g. 
[HHHI)) the lack of an underlying theory which explains 
the actual relevance of these corrections makes the ex- 
trapolation potentially uncontrolled. In particular, it is 
important to understand that one does not simply need 
to reach small densities, but it is also necessary to include 
in the fit only data in a region close to n = where a 
proper power-law behavior is present. Our investigation 
shows that this region can be very small, and densities 
of the order ot n ~ 0.1 are definitely outside and should 
not be included in the fit. 

For example, the different values of the critical temper- 
ature at unitarity calculated in Ref. |l3| and Ref. |l(| are 
likely explained by the lower minimal density used in Ref. 
[Hi , n m i n ~ 0.01, five times smaller than n m i n w 0.05 in 
Ref. [l(|, even if a different size extrapolation has been 
invoked to justify the discrepancy. 

We also notice that the linear fit to extrapolate T c at 
unitarity in Ref. [l3| includes data at density as high 
as n m 0.3, which are clearly not in the asymptotic 
regime according to the results in this paper and in Ref. 
(la ]. Also this effect is likely to affect the extrapolation 
and, as a general remark, only the data at the lowest 
should be used in the fit procedures in order to extract 
precise estimates for the dilute Fermi gas free from un- 
controlled errors. 



VI. CONCLUSIONS 

In this paper we studied the connection between the 
BCS-BEC crossover in dilute gases and in lattice models 
at finite density. More specifically we investigated within 
mean-field how several relevant observables in different 
lattice models all converge to those of a dilute gas as the 
density is reduced along suitable trajectories with con- 
stant gas parameter kpa s . Comparing with the solution 
of a dilute Fermi gas in continuum space, we analyzed 
quantitatively how the properties of lattice systems con- 
verge to the dilute limit. 

In analogy with the finite-size effects, which can be 
minimized through the knowledge of the expected de- 
pendence on the size of the system, a theoretical under- 
standing of the power-law density corrections to the di- 
lute regime is a crucial step towards a proper zero-density 



extrapolation. The mean-field results presented in this 
paper, though approximate, are meant to be a guideline 
to control the extrapolation. Moreover, mean-field re- 
sults are free from the interplay with finite-size effects, 
since our approach directly deal with the problem in the 
thermodynamic limit. 

Our results clearly show that, despite the universal na- 
ture of the dilute limit, finite density corrections depend 
on the specific choice of lattice. Moreover, these correc- 
tions persist down to extremely low densities, making any 
extrapolation based on lattice model at finite density very 
hard unless data at extremely low-densities are available. 
The leading order corrections to all the relevant observ- 
ables at small densities follow a n 1 / 3 power-law behavior 
for every lattice studied here and in Ref. [l8j |. The pref- 
actors instead depend on the specific lattice and vary 
over several order of magnitude as the lattice parameter 
spans the BCS-BEC crossover, being larger on the BEC 
side. The n 1 / 3 behavior also holds for a model with the 
same DOS of free fermions except for a sharp energy-cut- 
off which enforces the reduced translational symmetry of 
the lattice. In this case we also derived an exact ana- 
lytical expression for the leading-order prefactors in the 
whole crossover. 

This result shows unambiguously that the presence of a 
high-energy cut-off A associated with the broken Galilean 
invariancc affects the properties of the system down to 
arbitrarily low finite densities. This finding contradicts 
the common belief -valid only in the weak-coupling BCS 
regime- that only the low-energy behavior up to the 
Fermi energy Ep cx ?i 2 / 3 <C A determines the leading- 
order correction to the dilute limit. Comparing different 
lattices we show that also the differences of the DOS at 
finite energy far from the Fermi level contribute to the 
leading-order corrections. As a consequence, the prefac- 
tors of the leading density corrections are non universal 
and there is no simple recipe to cancel or even minimize 
them. 

The relevance of the DOS contrasts the idea that the 
leading-order corrections depend only on the low-energy 
scattering properties. This rules out the possibility to 
cancel the leading-order density corrections tuning the 
effective range to zero as proposed in [l7| . and seems to 
exclude a simple way to express the corrections in terms 
of a single parameter. 

The inclusions of fluctuations beyond mean-field is 
crucial for a quantitative description of the BCS-BEC 
crossover both in the dilute gas and in lattice models. 
However, fluctuations can even increase the sensitivity 
of the system to density variations, as found in Ref. [18[ 
comparing static and dynamical mean-field theory for lat- 
tice models at unitarity and in [30| where the relevance 
of density corrections to extrapolate the critical tempera- 
ture is investigated including Gaussian fluctuations. The 
same enhancement of lattice-dependent features is ex- 
pected to happen in finite-temperature data with respect 

to t = o Hi Hoi. 

During completion of the present work, we became 
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aware of a very recent Monte Carlo study|32j which shows 
a sizable effect of reducing the fermionic density from 
n ~ 0.1 to n ~ 0.04, in line with our results. 
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